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Based on the generalized valence bond picture, a Schwinger boson mean field theory is applied 
to the symmetric SU(4) spin-orbital systems. For a two-leg SU(4) ladder, the ground state is a 
spin-orbital liquid with a finite energy gap, in good agreement with recent numerical calculations. 
In two-dimensional square and triangle lattices, the SU(4) Schwinger bosons condense at (n/2, 7r/2) 
and (-71-/3,71-/3), respectively. Spin, orbital, and coupled spin-orbital static susceptibilities become 
singular at the wave vectors, twice of which the bose condensation arises at. It is also demonstrated 
that there are spin, orbital, and coupled spin-orbital long-range orderings in the ground state. 
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Recently the properties of spin systems with orbital de- 
generacy are attracting a lot of attentions .El Several spin- 
orbital models are proposed in various kinds of materials, 
such as C 60 ,H NiLi0 3 J Na 2 Ti 2 Sb 2 oJ4 and LaMn0 3 Ja In- 
terplay of spin and orbital degrees of freedom produces 
not only new magnetic structure phases but also novel 
quantum ordered and disordered states such as orbital 
density wave and spin-orbital liquids n |-A simplified and 
symmetric model for these systems isElQ 

H = Js ( 2S * ' S <+« + V2) (2T< • T i+S + 1/2) , 

1 i,S 

(1) 

where the operators S and T are SU(2) Pauli matrices 
for the spin and orbital degrees of freedom, respectively. 
The vector S points to the nearest neighboring sites. It is 
already known that the model possesses SU (4) symmetry, 
and can be derived from a quarter-filled electronic model 
with two-fold orbital degeneracy by ignoring the Hund's 
rule couplings and taking the large on-site Coulomb inter- 
action. High symmetry in this spin-orbital model means 
strong correlations between spins and orbitals. In an 
SU(4) symmetric state, the correlation functions for fif- 
teen generators of the SU(4) Lie group are isotropic. The 
spin, orbital, and coupled spin-orbital degrees of free- 
dom must be treated on equal footing.tffl Over last a 
few years this model was studied extensively. In one 
dimension (ir^_iOs_fairly understood analytically and 
numericallyJjoOE3£3 The ID model can be solved by 
means of Bethe ansatz, and its ground state is described 
by a gapless spin liquid, similar to the SU(2) Heisenberg 
model. In two dimension (2D), it is relatively less un- 
derstood. Li, et alB first argued that an SU(4) singlet 
plaquette state contains at least four sites, and a col- 
lection of such SU(4) singlets may lead to a spin liquid 
state. Thia-jpicture is realized very weiLin two- leg lad- 
der model,li-3 and some solvable modelsJla So far it is not 
clear whether such an idea can be realized ip-2D. Except 
numerical diagonalization on small clusterst£l and series 
expansion,LJ there has been no solid results as whether 
the ground state is long-range ordered or a spin-orbital 
liquid. 

In this paper, the properties of the generalized valence 



bond state consisting of the SU(4) singlets are discussed 
in detail for the model (Eq.(|l|)) in a SU(4) Schwinger 
boson mean field theory. The ground state of the sys- 
tem can be either ordered or disordered, depending on 
the dimensionality and lattice topological structure. In a 
two-leg ladder system, we find that the ground state is a 
spin liquid state with a finite energy gap, which is in good 
agreement with recent numerical calculations by van den 
Bossche.EJ In square and triangle lattices, the Schwinger 
bosons condense at zero temperatures, i.e., the Bosc- 
Einstein condensation (BEC) occurs and is identified as 
the indication of LROs in the ground state. The relation 
between BEC and LROs is illustrated explicitly in spin, 
orbital, and spin-orbital static transverse susceptibilities, 
which become singular at the wave vectors (7r/2, 7t/2) and 
(tt/3, 7t/3) for the respective square and triangle lattices, 
leading to finite staggered magnetizations for the spin, 
orbital, and coupled spin-orbital densities in the thermo- 
dynamic limit. Thus, three Goldstone modes are found. 
According to the calculated susceptibilities, the spin, or- 
bital, and spin-orbital LROs may coexist, but the direc- 
tion of spontaneous symmetry breaking will determine 
the properties of the ground state. 

In general, for the spin- 1/2 system with double or- 
bital degeneracy, there are four local states on each site 
i according to the eigenvalues of Sf and r|: |1) = 
1 + 1/2, +1/2), |2) = |-l/2, +1/2), |3) = |+l/2,-l/2), 
|4) = | — 1/2,— 1/2). Four Schwinger bosons can be in- 
troduced to describe these four states: |/x) = aj, |0) 
where |0) is the vacuum states and \i = 1,2,3,4. There 
has to be imposed a local constraint, X^=i a !^ a v — 
1 on each lattice site. The permutation operator 
P i:j = (2S, • Sj + |) (2Ti • T 3 + |) is to exchange the 
two states on the sites i and j, Pij \ifi,jv) = \w,jy). 
Moreover, can be expressed in terms of the four 

hard-core bosons as P io = J2 fJ ., u a l l j,^ a ]iy a 3tJ-- An SU(4) 
singlet is defined by SU^iJ^kJ) = Y^,^,u,y,s^^,j,S 

a \^ a jv a \~j a \6 |0)> wnere ^n,v.~).s is an antisymmetric ten- 
sor. When the model Hamiltonian Eq.(Q) has only four 
lattice sites, an SU(4) singlet is always the lowest energy 
state for Js > 0. According to the group theory, the 
SU(4) symmetric state for a lattice with An sites (n is 
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integer) can be regarded as a linear-combination of all 
states consisting of n SU(4) singlets.Ej This is a general- 
ization of the Anderson's resonating valence bond (VB) 
stateEa from the spin SU(2) system to the SU(4) system. 
As is well-known in the Heisenberg model, a short-range 
VB state may describe a spin liquid state with a finite 
energy gap,E£l and a long-range VB state may possess 
antiferromagnetic LROHj A Schwinger boson mean field 
theory based on the short-range VB state was proposed, 
for the spin SU(2) systems by Auerbach and Arovas,Eil 
which successfully describes either ordered or disordered 
quantum states. Very recently, the theory was applied 
to the spin-orbital systems by the present author and his 
collaborator^ 

To realize the generalized VB state in the SU(4) spin- 
orbital system, the model Hamiltonian is rewritten as 

i, 5,^,1/ 

+E A *(E a L a ^- 1 ) + 2 i v A E J <5 ' (2) 

i fj,— 1 8 

where Aij- tfljU = a^aju — a^a^ and N\ is the total 
number of lattice sites. Antisymmetric operators j ;MiJ/ 
are introduced for the purpose of the mean field cal- 
culations. The following theory is limited to the case 
J$ > 0. The local Lagrangian multiplier is also used to 
impose the local constraint for the hard-core bosons on 
average. In the mean field approximation we will take 
all Ai = A. The thermodynamic averages of the op- 
erators Aij^ u are defined as the VB order parameters, 
{Aij+B;ii,v) = — 2«A jU ,„(5), which are odd functions with 
respect to either the indices fi, v or the vector direction 
5. In the momentum space, we define 

7^(k) = 2*£j a A /H ,(tf)ex P (ik.«). 
s 

We define an eight-component spinor 

*k = (°kl' a k2' a k3> a k4' a -kl, a_ k 2, Ct-k3, a-k4)- 

By utilizing the Pauli matrices a a (a = x, y, z), the de- 
coupled mean field Hamiltonian is thus written in a com- 
pact matrix form, 

H = ^E $ k( A -^® B ( k )) $ k+f (3) 

k 

where £ /N A = J2 S J a A%,(S) -3X+^Es ft; 

( 7l2 (k) 7l3 (k) 7i 4 (k)\ 
Rn ^ _ -7i 2 (k) 723 (k) 7 24 (k) 

B w - _ 7i3(k) _ 723(k) o 734(k) , 

\-7u(k) -7 24 (k) -7 34 (k) / 

where N\ is the number of lattice sites. Considering 
the symmetry in the HamiltonianJ23 there exists two sets 



of solutions: (I). 7l2 (k) = 7 34 (k), 7l3 (k) = -7 24 (k), 
7l4 (k) = 723 (k); (II). 7l2 (k) = -7 34 (k), 7i 3 (k) = 
7 24 (k), 7 14 (k) = — 723 (k). The physical reason is that 
an SU(4) singlet plaquette state SU4,(i,j,k,l) contains 
four creation operators with different sites and indices, 
and A]- ,,,, contains only two sites and two indices. To 
form such an SU(4) singlet plaquette state, the four cre- 
ation operators on different sites should have different 
indices. The relations in the solution (I) and (II) reflect 
these properties. In each SU(4) singlet plaquette it con- 
sists of two configurations. Each configuration consists 
of two spin singets and two orbital singlets. The two 
configurations are degenerated, but not orthogonal. The 
double degeneracy of the solutions may be related to the 
properties of SU(4) singlet plaquettes. So these relations 
will help us to construct the wave function of the gen- 
eralized SU(4) VB state. We first focus on the solution 
(I), and will discuss the results of the solution (II). The 
single-particle Green function is given by 

G(k, iujn) = [iuj n VL\ — A + iu y ® B(k)] -1 
iunflx + A + 7 12 (k)» 2 + 7 13 (k)r? 3 + 7 13 (k)^ 4 

(iw n )2 - W(k)2 U 

where = <r z ® Co ® <tq, fi 2 = <x y ® cto ® cr y , = 
g v (8 <Jy ® 0>, an d f^ 4 = <J y <E> <J y <E> <J x ■ There is only one 
four fold degenerate quasiparticle spectrum: 

"00 = \j>?- [ 7l2 (k)+7 13 (k)+7 14 (k)], (5) 
from which the free energy for the system is evaluated 

F = \ E ln [l " exp[-/3^(k)] + 2 "(k) + So. (6) 

P k k 

The saddle point equations are thus derived by minimiz- 
ing the free energy with respect to the mean field vari- 
ables A and A M „(<5). We can deduce the VB order pa- 
rameters according to the symmetry of A flL ,(S) and of 
lattice. 

Now we apply the general formalism to the SU(4) 
model on several lattices. We first study the two-leg lad- 
der model. Recent numerical study has shown that its 
ground state is a spin-orbital liquid with a finite energy 
gapjifl which can be regarded as a realization of SU(4) 
plaquette state, or a short-range generalized VB state. 
Here the isotropic case is considered J|| = J± = J so that 
we introduce two sets of VB order parameters: A^ v (x) 
along the ladder and A Atl/ (y) along the rungs. Special at- 
tention should be paid for the direction along the rung. 
The momentum along the rung has two discrete values. 
The spectra are given by 

w±(k) = yj\ 2 - I6J 2 A 2 (sink±7/) 2 (7) 

where Ajj = A? 2 (x) + Aj 3 (x) + A 2 4 (x) and 77 = 
A^ J /(y)/2A /JI/ (a;) is determined by the mean field 
equationscJ. Two branches of spectra have a relation: 
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.(k) 



-k). The saddle point equations are deter- 



mined by minimizing the free energy with respect to the 
mean field variables 



dk A 

2^w+(k) 
dk (sink + 77) 2 
2?r w+(k) 
dk sin k + rj 
2tt w+(k) 



[2n fl (w + (k))+l] 
[2n fl (w+(k))+l] 
[2n B (w+(k)) + l] 



3 
2' 

1 + 2?7 : 



(8a) 



2 J 



J' 



{8b) 
(8c) 



where n B (x) is the distribution function for bosons: 
n B (x) = l/[ex$\px] - 1] and /3 = l/k B T. At T = 0, 
ns(w+(k)) = when w+(k) ^ 0. The numerical cal- 
culations give rise to A = 2.26574J, A|| = 0.35339, 
and rj = 0.55146. The minima of spectra is at k* = 
7r/2 for uj + (k), and at —ir/2 for w_(k): min(o;±(k)) 

= wA 2 — 16J 2 A 2 (1 + rf) 2 . From the dynamic suscep- 



tibilities of spin, orbital, and spin-orbital operators, 
we find that there is a finite energy gap: Aqap — 
2min(w±(k)) = 1.138 J. For the second set of solution 
(II), it also produces the same numerical results, degen- 
erated with the first set of solutions, van den Bossche et. 
alS-i studied this SU(4) ladder model up to 16 sites by an 
exact diagonalization method, and a finite energy gap has 
been found in a singlet-multiplet excitation: A = 1.09 J. 
The two values are in excellent agreement. A local min- 
imum at k = 7r/2 in the quasiparticle dispersion was 
also observed. Moreover, the ground state has also been 
found to have a two-fold degeneracy in the thermody- 
namic limit, consistent with the two solutions of ours. 
All these facts can be regarded as strong support for 
our present theory. It is worth noting that the theory 
may fail for one-dimensional chain to predict a small en- 
ergy gap due to the ignorance of the topological terms 
or "Pontryagin index" , which can destroy the Haldane's 
gap. The same problem was encountered in the SU(2) 
theory, and was discussed extensively in one-dimensional 
spin 1/2 systems. Ell 

Next we come to study 2D lattices. Let us consider a 
square lattice first. In this case, we still assume isotropic 
couplings J x — J y . The spectrum for the bosonic quasi- 



particles can be written as 



w(k) = J A - 16 J 2 A 2 (sinks + sink,,) 2 



(9) 



with A 2 = A 2 2 + A 2 3 + A 2 4 . The minimum of the energy 
spectra occurs at k* = (tt/2,tt/2). The saddle point 
equations are given by, 

c^k \ 3,, 
(2^^(k) ^{k) ) + l] = -(10) 



dk (sin k x + sin k y ) z 



[2n B ( W (k)) + l] = ^,1) 



(2tt) 2 w(k) 

where the coordinate number for square lattice is z — 4. 
In the present theory the number of bosons in the di- 
agonalized Hamiltonian is not equal to the number of 



the hard-core bosons since the Bogoliubov transforma- 
tion changes the number of bosons. The quasi-particle 
number is determined by solutions of the saddle point 
equations self-consistently. At T = the bosons may 
condense, i.e., the BEC occurs. The minima of the 
boson spectra a>(k) are at k = ±k*. Since the dis- 
tribution function becomes singular at w(k*) = in 
Eq. (|Io[) , we have to introduce a finite quantity 60 = 
2[n B (uj(k*)) + n B (uj(-k*))} / [N A uj(k*)} such that the 
saddle point equations have physical solutions. The non- 
zero 60 m ay be used as the order parameters for long- 
range correlations for the spin, orbital, and coupled spin- 
orbital densities. The saddle point equations are solved 
numerically. When T > 0, 60 = 0. At T = 0, we have 
A = 8JA = 3.05905J and b = 0.1068. The non-zero b 
indicates that the BEC occurs on a square lattice. 

Another typical 2D lattice is the triangle lattice. For 
example, NiLiOa has a 2D triangle lattice-structure, and 
was modeled as a spin-orbital systemO'EEl. Usually the 
quantum frustration is anticipated to make quantum fluc- 
tuations more stronger. Topologically, we can distort a 
triangle lattice into a square one by introducing a finite 
diagonal coupling J x +y, and take another diagonal cou- 
pling J x - y = 0. In the present theory we have three sets 
of order parameters: A M „(x), A M „(y), and A M „(x + y). 
Focus on the isotropic case J x = J y = J x + y , then the VB 
order parameters have a relation: A M „(x) = A M „(y) 



A M1/ (x 



= A 



//.// • 



The spectra for the bosonic quasi- 



particles is thus given by 



w(k) = y\ - 16 J 2 A 2 (sinks +sink y + sin(k x + k y )) 2 , 

(12) 

with A 2 = A^ 2 + A^3 + A 2 4 . The minimum of the spectra 
occurs at k* = ±(7r/3, 7r/3), which can be shifted away 
if the couplings are anisotropic. A similar set of saddle 
point equations are obtained by minimizing the free en- 
ergy in Eq. (g). Numerically solving the self-consistent 
equation at T = gives rise to A = 6\/3</A = 3.57878 J 
and b Q = 0.155345. Since w(k*) = 0, the BEC also ap- 
pears on an isotropic triangle lattice. The coupling J x + y 
does not enhance the quantum frustration to suppress the 
BEC completely. The role of J x + y is to force the minimal 
point from ±(7r/2,7r/2) for J x + y — to ±(7r/3,7r/3) for 
Jx+y — Jx = Jy- For a more detailed calculations, it is 
found that the wave vector k* changes continuously as a 
function of J x +y/Jx- 

In the usual Schwinger boson mean field theory the 
BEC is identified as the long-range correlations between 
the SU(2) spins. To establish the relation between the 
BEC and the long-range orders in the present SIK4) 
ground state, we calculated the static susceptibilitiesEj 



xx(q) 



— V 

Ma ^ 



16A A 



xTr{r2 x G(k + q,T = Q-)n x G{k,T = 0+} , (13) 

where X = S for spin Sf, X = T for orbital T?, and 
X = ST for the coupled spin-orbital 2S|Tf. These 
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three operators can be expressed in terms of spinors, 
$Tf2x$i/4, where tl s = °o ® &o ® , = oo® *® ^) 
and f^5T = Co ® o"« ® °z ■ From the single particle Green 
function Eq.(0), we can calculate the static susceptibili- 
ties, 



When A 2 2 



A 2 

^13 



A 2 

^14 



Xs(Q)/Na 
Xt(Q)/Na 
Xst(Q)/Na 



& 2 7 ? 3 (q = 2k*), 
& 2 7 2 2 (q = 2k*), 
& 2 7 2 4 (q = 2k*), 



(14a) 
(14b) 
(14c) 



which become singular when q = Q = 2k* and higher 
order terms can be ignored. The correlation functions 
are proportional to the number of lattice sites once 
bo =/= 0. These properties are characterized by long- 
range correlations at the wave vector Q. In the ther- 
modynamic limit, the corresponding magnetizations be- 
come mx = \J Xx(Q)/Nk, which depend on the values 
of VB order parameters, A 2 2 , A 2 3 , and A 2 4 , respectively. 
The long-range order is thus at Q = (jr, ir) for a square 
lattice ((7r, 7r) and (— n, — n) are the same vector), and 
Q = (27r/3,27r/3) or (-2^/3,-2^/3) for a triangle lat- 
tice. The two vectors ± k* correspond to one state for 
a square lattice, but two equivalent states for a triangle 
lattice. From the second set of mean field solutions, we 
find that the above relations remain if we make a permu- 
tation between S and T. Therefore these two solutions 
(I) and (II) are degenerated. We believe that the double 
degeneracy of the ground states observed in our theory 
is not a result of the mean field theory, and may have 
a deep physical origin. The singularity in static suscep- 
tibilities also reflects the fact that the collective modes 
are gapless Goldstone modes. The SU(4) system may 
have at most three Goldstone mode when the symmetry 
is broken spontaneously. 



A 2 /3, we have \s — Xt — 
Xst- These relations are in agreement with the isotropic 
correlations for spin, orbital and coupled spin-orbital in 
our SU(4) invariant ground state. In such a mean field 
theory, the operator S z , T z , and 2S Z T Z play an equal 
role. When spin and orbital density operators have long- 
range correlations, and the coupled spin-orbital density 
operators also have long-range correlations with the same 
wave vectors. An interesting observation is that t 
ground state energy depends only on the parameter A 2 . 
It may contain some new states, which is determined by 
the direction of spontaneous symmetry breaking.c3 Actu- 
ally, we have the freedom to choose the direction of the 
spontaneous symmetry breaking in the thermodynamic 
limit. For example, an infinitesimal external staggered 
magnetic field along the spin z-direction may induce a 
ground state with A 2 2 = A 2 and A 2 3 = A 2 4 = 0, where 
a magnetic long-range order appears only. An infinites- 
imal John- Teller distortion may induce a ground state 
with A 2 3 = A 2 and A 2 2 = A 2 4 = 0, where an orbital 
long-range order shows up. It is also possible for two or 
more types of long-range orders to coexist in the a single 
ground state. 

In conclusion, based on a generalized valence bond 
state picture, a Schwinger boson mean field theory is de- 
veloped for the symmetric SU(4) spin-orbital systems, 
showing that the ground state for a two-leg ladder model 
is a spin-orbital liquid with a finite energy gap in low 
energy excitations, and the ground states for square and 
triangle lattices possess spin, orbital, and coupled spin- 
orbital long-range orderings. 
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no SU(4) symmetry. To keep the SU(4) symmetry of the 
state and the isotropic correlation functions we have to 
make some restriction on the 6 order parameters. We take 
w+(fc) = ui-(k) = ui(k) such that the the spectra become 
four-fold degenerated. The equality gives the two sets of 
the solutions. 

r\ is independent of [iv which can be derived from the mean 
field equations if we list all equations for A p „(<5). 
The spectra in Eqs. (9) and (11) have no rotation 
symmtries of square and triangle lattices. Under the rota- 
tion symmetry we have several other solutions which give 



the same physical consequences. It is worthy noting that 
the dispersion relations for Schwinger bosons are not ob- 
servables and physical observables are related to the spin 
operators explicitly. 
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